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Abstract.  This  paper  presents  an  analysis  on  the  space/time  statistical  thermal  structure 
in  the  Yellow  Sea  from  the  Navy's  Master  Observation  Oceanography  Data  Set  during 
1929-1991.  This  analysis  is  for  the  establishment  of  an  Optimum  Thermal  Interpolation 
System  of  the  Yellow  Sea  (a  shallow  sea),  for  the  assimilation  of  observational  data  into 
coastal  <r  coordinate  ocean  prediction  models  (e.g.,  the  Princeton  Ocean  Model),  and  for 
the  design  of  an  optimum  observational  network.  After  quality  control  the  data  set 
consists  of  35,658  profiles.  Sea  surface  temperatures  at  50%  and  80%  water  depths  are 
presented  here  as  representing  the  thermal  structure  of  surface,  middepth,  and  near- 
bottom  layers.  In  the  Yellow  Sea  shelf  the  temporal  and  spatial  signals  fluctuate  according 
to  the  Asian  monsoon.  Variation  of  surface  forcing  from  winter  to  summer  monsoon 
season  causes  the  change  of  the  thermal  structure,  including  the  decorrelation  scales.  Our 
computation  Shows  that  the  seasonal  variation  of  the  surface  horizontal  decorrelation  scale 
is  around  90  km  from  158  km  in  w  inter  to  251  km  in  summer  and  the  seasonal  variation 
of  the  surface  temporal  decorrelation  scale  is  around  2.4  days  from  14.7  days  in  winter  to 
12.3  days  iri  summer.  The  temporal  decorrelation  scale  increases  with  depth  in  both 
summer  (evident)  and  winter  (slight).  The  near-bottom  water  (tr  -  0.8)  has  the  longest 
temporal  scale  in  summer,  which  could  be  directly  related  to  the  existence  of  the  Yellow 
Sea  Cold  Water  throughout  the  summer  in  the  middle  of  the  Yellow  Sea.  The  temporal 
and  spatial  decorrelation  scales  obtained  in  this  study  are  useful  for  running  optimum 
interpolation  models  and  for  designing  an  optimum  observational  network.  The  minimum 
sampling  density  required  to  delect  thermal  variability  in  the  Yellow  Sea  shelf  would  be 
50—80  km  and  4-6  day  intervals  per  temperature  measurement  with  the  knowledge  that 
the  Subsurface  features  will  also  be  adequately  sampled. 


1.  Introduction 

Various  ocean  systems  such  as  fronts,  eddies,  and  water 
masses  have  different  temporal  and  spatial  scales.  These  scales 
feature  a  system's  life  span  and  spatial  extent  both  horizontally 
and  vertically.  Horizontal  scales  for  temperature  fields  in  deep 
ocean  basins  have  been  extensively  studied:  for  example.  White 
ct  al.  [1982J  identified  spatial  correlation  scales  in  the  western 
Pacific  of  about  60(1  km  in  the  tropics  (south  of  17.5  N  )  and  300 
km  in  the  subtropics  (north  of  1 7.5'N),  while  Ozsoy  et  al.  (  1  989] 
found  the  spatial  scales  to  be  200-250  km  in  the  Levantine 
Basin  of  the  Mediterranean  Sea. 

Temporal  and  horizontal  scales  in  the  continental  shelf,  and 
especially  in  the  region  of  freshwater  influence,  have  been  levs 
studied.  One  major  reason  is  that  due  to  great  depths  the  deep 
ocean  temperature  profiles  tend  to  follow  the  pattern  of  mixed 
layer,  a  thermodine.  and  a  deeper  layer  slowly  decreasing  in 
temperature  w'ith  depth.  Shallow  water  docs  not  consistently 
mimic  this  pattern.  It  may  range  from  the  classical  profile  to 
completely  isothermal.  The  coastal  water  is  largely  affected  by 
bathymetry,  river  runoff,  internal  waves,  and  tides.  Therefore 
the  temporal  and  horizontal  scales  are  also  under  the  strong 
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influence  of  these  forcing  factors  and  should  have  different 
vertical  structure  from  the  deep  water  [  Wells.  1994J.  For  deep 
waters  the  vertical  variations  of  the  temporal  and  horizontal 
scales  are  usually  obtained  from  the  temperature  fields  at  dif¬ 
ferent  depths  (e.g..  100,  200.  and  300  m  from  White  el  al. 
1 1982])  or  from  the  sea  surface  temperature  (SST)  field  and 
the  depth  of  certain  temperature  (e.g..  depth  of  20°C  from 
Spriniall and  .Mom  [  199|  |).  Neither  scheme  is  practical  for  the 
Yellow  Sea  region  because  of  a  single-layer  structure  appear¬ 
ing  in  winter,  which  makes  the  depth  of  certain  isotherm  non- 
representative  for  the  subsurface  thermal  fields,  and  because 
of  a  strong  influence  of  horizontally  vary  ing  bathymetry  on  the 
shallow  water,  which  causes  the  temperature  fields  at  different 
depths  to  be  nonrepresentative  of  subsurface  fields.  For  deep 
water  these  depths  are  easily  determined  (Figure  la).  How¬ 
ever.  for  shallow  water  with  a  sloping  bottom,  any  chosen 
depths  will  lead  to  some  areas  (hatched  area  in  Figure  lb)  not 
covered  by  any  submodels.  Therefore  no  climatology  can  be 
constructed  for  these  areas.  In  order  to  solve  this  problem,  we 
use  tr  surfaces 

<T  =  -2/H  (1) 

to  build  up  climatology  and  statistical  structure.  Here  H  is  the 
water  depth,  and  z  is  the  vertical  coordinate.  Midwater  prop¬ 
erties  are  represented  by  a  =  0.5.  and  the  near-bottom  water 
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Figure  1.  Depths  for  building  temperature  climatology:  (a) 
horizontal  levels  for  deep  water,  (b)  horizontal  levels  for  con¬ 
tinental  shelf,  and  (e)  bottom  following  levels. 


features  are  portrayed  by  a  =  0.8  (Figure  lc).  The  analysis  on 
the  a  coordinate  system  is  a  benefit  for  initializing  popular  <r 
coordinate  coastal  models,  such  as  the  Princeton  Ocean  Model 
[Blum berg  and  McUor.  1987].  However,  the  analysis  on  the  it 
coordinate  system  has  disadvantage  during  weakly  forced  and 
stratified  periods.  The  deeper  n  levels  will  be  within  the  mixed 
layer  environment  in  shoal  areas  and  below  the  mixed  layer  in 
the  deeper  portions  of  the  basin. 

'Hie  U.S.  Navy's  Master  Oceanographic  Observation  Data 
Set  I  MOODS)  has  approximately  35.658  profiles  (during 
1929-1991 )  for  the  Yellow  Sea  shelf.  After  climatology  has 
been  established  from  these  profiles,  temperature  anomalies 
are  computed  by  subtracting  the  climatology  from  the  obser¬ 
vations.  An  autocorrelation  function  (ACF)  is  then  defined  by 


Tjtl  )  = 


7"(/0)7'(/0+  l)dla 


(-) 


where  7”  is  the  temperature  anomaly,  l„  denotes  the  indepen¬ 
dent  space, ’time  vectors  defining  the  location  of  points  in  a 
sampling  space  /,.  and  /  is  the  space  time  lag  and  v2  is  the 
variance.  Here  7|  is  computed  by  paring  the  anomalies  into  bins 
depending  upon  their  separation  in  space  time.  /.  The  values  of 
77  will  be  obtained  from  calculating  the  correlation  coefficient 
for  all  the  anomaly  pairs  in  each  bin  which  will  be  constructed 
for  Lite  combination  of  different  lags. 


2.  Characteristics  of  the  Yellow  Sea 

The  Yellow  Sea  is  the  middle  section  of  a  continental  sea 
system  (Figure  2).  To  the  north  is  the  Bohai  Sea  and  to  the 
\outh  are  the  East  China  Sea  and  the  western  Pacific  beyond. 
It  lies  between  China  and  the  Korea  Peninsula  and  covers  an 
area  of  approximately  295,000  km2.  There  arc  four  major  rivers 
that  flow  into  it.  The  Yangtze  River  is  located  to  the  southwest, 
the  Yellow  River  and  the  Liao  River  to  the  north,  and  the  Hun 
River  to  the  east.  While  the  Yetlow  Sea  covers  a  relatively 
large  area,  it  is  uniformly  shallow.  The  water  depth  over  most 
of  the  area  is  less  than  50  m.  The  deepest  water  is  confined  to 
a  north-south  oriented  trench  paralleling  the  Korean  coastline 


in  the  central  portion  with  a  maximum  depth  of  90  m.  The 
bottom  slopes  from  20  m  near  the  China  coast  to  its  southern 
boundary',  the  shelf  break,  at  150  m.  Within  50  km  of  the 
Korean  coastline  the  average  water  depth  is  20  m. 

The  Yellow  Sea  lies  within  the  Asian  monsoon  circulation 
area  and  experiences  hot.  humid  summers  and  cold,  dry  win¬ 
ters.  In  winter  the  Yellow  Sea  is  dominated  by  the  cold  north¬ 
erly  wind  associated  with  the  Siberian  high  (Figure  3a).  The 
wind  speed  reaches  a  maximum  in  February  with  a  value  of  35 
ms  near  Ihe  central  portion  and  of  28  m/s  over  the  whole 
Yellow  Sea  [Elms.  1990].  The  northerly  wind  bursts  create  a 
southward  sea  level  gradient  that  forces  bottom  water  to  flow 
northward  in  spurts  along  the  Yellow  Sea  trough  [Hsiii’h. 
1988].  The  January  surface  air  temperature  (SAT)  varies  from 
O'  to  8T  in  the  Yellow  Sea  (Figure  4a).  roughly  2r-6  C  cooler 
than  SST.  The  Yellow  Sea  surface  loses  heat  to  the  atmo¬ 
sphere.  The  upward  buoyancy  flux  at  the  air-ocean  interface 
(thermal  forcing),  together  with  the  strong  wind  stress  (me¬ 
chanical  forcing),  generates  turbulence  and  mixes  the  surface 
water  with  the  deeper  water.  The  mixed  layer  is  at  its  deepest 
(it  usually  fills  the  whole  water  column)  during  winter  owing  to 
both  convection  and  wind  mixing  by  the  strong  northeast  mon¬ 
soon  winds.  The  change  of  seasons  begin  in  March  when  the 
surface  air  temperatures  are  5°C  warmer  than  in  February. 
Rapid  weakening  of  the  Siberian  high  progresses  into  April.  In 
late  April  the  atmospheric  polar  front  has  moved  northward 
toward  Korea  with  warm,  moist  air  following  behind.  Numer¬ 
ous  front-driven  events  occur,  making  late  April  and  May 


Figure  2.  Bottom  topography  of  the  Yellow  Sea  and  the  sur¬ 
rounding  regions.  The  data  are  obtained  from  the  U.S.  Naval 
Oceanographic  Office  DBDB5  world  bathymetry  database. 
Depths  are  in  meters. 
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highly  variable  in  terms  of  winds,  clouds,  and  precipitations. 
The  Yellow  Sea  surface  quickly  transits  from  winter  to  spring 
temperatures  with  an  average  increase  of  ICC. 

In  late  May  and  early  June  the  atmospheric  low-pressure 
system  begins  to  form  in  the  north  of  the  Yellow  Sea  and  to 
migrate  toward  the  west  over  Manchuria  in  late  June,  setting 
up  the  southwest  monsoon  that  dominates  the  summer 
months.  The  Manchuria  Low.  associated  with  the  atmospheric 
high-pressure  system  in  the  southeast  (the  Bonin  Highl.  pro¬ 
duces  southerly  winds  carrying  warm,  moist  air  over  the  Yellow 
Sea  (Figure  3b).  The  summer  S  AT  is  quite  uniform  and  around 
24  -26  C  (Figure  4b).  It  is  usually  1.5  -2  C  warmer  than  SST 
[ Van  Loon,  1984  J.  The  warm  air  and  the  strong  downward  net 
radiation  cause  downward  heat  flux  and  stabilize  the  upper 
layer  of  the  water  and  cause  the  surface  mixed  layer  to  shoal. 
Associating  with  the  shallow  mixed  layer  a  closed  surface  cy¬ 
clonic  circulation  occurs  above  the  thcrmocline.  which  is  about 
13  m  in  August,  and  may  reach  a  speed  of  15  cm/s  \Bartz.  1972]. 
Below  the  thcrmocline.  there  is  a  cold  water  mass,  commonly 
referred  to  as  the  Yellow  Sea  Cold  Water  (YSCW)  mass,  that 


Figure  3.  Mean  atmospheric  surface  circulations  in  the  vi- 
cinitv  of  the  Yellow  Sea  for  (a)  February  and  (b)  June  |from 
Lunghill.  1976]. 


Figure  4.  Mean  surface  air  temperature  in  the  vicinity  of  the 
Yellow  Sea  for  fa)  January  and  (b)  July. 


remains  unchanged  and  nearly  motionless  throughout  the  sum¬ 
mer  ]/_/  utul  Yuan,  1992).  October  is  the  beginning  of  the 
transition  back  to  winter  conditions.  The  southerly  winds  have 
weakened,  letting  the  sea  surface  slope  reestablish  toward  the 
winter  pattern  again.  The  SST  steadily  decreases  from  October 
to  January. 

3.  The  MOODS  Data  Set 

The  temperature  profiles  (1929-1991)  used  for  this  study 
were  taken  from  MOODS.  There  are  approximately  50.00(1 
temperature  profiles  located  between  119-128  0  and  25  - 
40  N.  which  encompass  the  Yellow  Sea  and  the  northern  East 
China  Sea.  The  data  were  screened  to  restrict  the  data  set  to 
those  profiles  that  are  located  on  the  continental  shelf.  This 
reduces  the  data  set  to  35.658  profiles  ( Figure  5).  To  investi¬ 
gate  the  seasonal  variation  of  the  temporal  and  spatial  scales. 
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LongtuOe  (E) 

Figure  5.  Distribution  of  MOODS  profiles  in  this  stud}. 

the  database  was  binned  into  four  seasons.  The  seasons  were 
defined  according  to  the  convention  of  the  Naval  Oceano¬ 
graphic  Office  for  the  Yellow  Sea.  January.  February,  and 
March  constitute  winter:  April.  May.  and  June  constitute 
spring:  July,  August,  and  September  constitute  summer:  and 
October.  November,  and  December  constitute  fall.  Table  1 
shows  the  partition  of  the  data:  w  inter  has  the  least  profiles  and 
summer  has  the  most. 

Fite  main  limitation  of  the  MOODS  data  Is  their  irregular 
distribution  in  time  and  space.  Certain  periods  and  areas  arc 
over  sampled,  while  others  lack  enough  observations  to  gain 
any  meaningful  insights.  Vertical  resolution  and  data  quality 
are  also  highly  variable  depending  much  on  instrument  type 
and  sampling  expertise.  There  is  a  data  sparse  area:  the  eastern 
coastal  region  of  China  I  Figure  5  ).  The  period  of  1963-1977  is 
found  to  have  relatively  large  number  profiles,  averaging  315 
per  season.  An  example  of  temporally  uneven  distribution  can 
be  seen  from  Figure  (>.  which  indicates  a  number  of  tempera¬ 
ture  observations  during  January-June  1977.  There  tire  fre¬ 
quent  15-20  day  gaps  of  no  observations  in  the  whole  Yellow 
Sea.  Spatial  and  temporal  irregularities  along  with  the  lack  of 
data  in  certain  regions  must  be  carefully  weighted  in  order  to 
avoid  statistically  induced  variability.  Originally,  it  was  hoped 
that  there  would  be  sufficient  data  to  compute  horizontal  and 
temporal  decorrelation  settles  on  a  yearly  basis  and  to  deter¬ 
mine  a  trend.  .After  analyzing  the  location  and  histogram  of  the 
data,  it  was  obvious  that  statistically  significant  results  cannot 
be  obtained  on  a  yearly  basis.  Analysis  was  done  on  a  mean 
seasonal  basis  using  the  data  for  all  years. 


Table  1.  Number  of  Temperature  Profiles  in  the  Database 
lor  Each  Season 


Number  of 

Season 

Profile- 

Winter 

739] 

Spring 

9868 

Summer 

99 15 

Fall 

S4S4 

4.  Establishment  of  Seasonal  Climatology 

Building  climatological  fields  from  the  MOODS  database  is 
not  a  straightforward  task,  and  in  shallow  water  the  problems 
are  magnified.  Most  of  the  problems  are  related  to  lack  of  data. 
Without  enough  data,  individual  observations  force  artificial 
gradients  to  occur.  In  addition,  the  sloping  bottom  poses  dif¬ 
ficulties  with  interpolating  techniques.  If  the  temperature  val¬ 
ues  need  to  he  unbiased,  an  appropriate  background  field 
(seasonal  climatology)  must  be  subtracted.  Three  major  steps 
were  used  in  establishing  the  seasonal  temperature  climatol¬ 
ogy:  (1 )  binning  the  data  seasonally,  (2)  interpolating  the  tem¬ 
perature  profiles  into  the  desired  depth,  and  (3)  gridding  the 
temperature  with  the  Minimum  curvature  SPline  (M1SP)  In¬ 
terpolator  algorithm  developed  at  the  Naval  Oceanographic 
Office.  This  algorithm  is  similar  to  the  Cubic  Spline  Interpo¬ 
lator  (c.g..  as  given  by  in  MAT! AH  [1992]).  During  the  M1SP 
interpolation  we  apply  a  low-pass  filter  to  the  observed  tem¬ 
perature  values  at  each  <r  level  of  interest.  That  is.  the  inter¬ 
polated  background  temperature  at  a  grid  was  a  weighted 
linear  combination  of  the  adjacent  measured  values  at  the 
same  <r  level.  Tile  choice  of  the  decay  scale  for  the  Gaussian 
weights  of  the  filter  was  somewhat  arbitrary  except  that  we 
wanted  it  to  be  much  longer  than  the  decorrelation  scale  of 
ACT  We  used  a  decay  scale  of  5  for  the  filler  weights  which 
gave  a  very  smooth  background  field. 

The  mean  temperature  field  at  tr  =  U  (Figure  7)  is  more 
uniform  horizontally  in  summer  than  in  winter.  The  mean 
temperature  fields  at  tr  -  0.5  t  Figure  8)  and  it  =  0.8  ( Figure  9) 
mimic  each  other  from  season  to  season.  The  only  difference  is 
that  the  temperature  reduces  near  2°C  in  summer  from  a  -  0.5 
to  it  =  0.8.  representing  12-24  m  in  depth  change  (Figure  2). 
The  winter  subsurface  (<r  =  0.5.  0.8)  mean  temperature  fields 
closely  mirror  the  mean  SST  field  (single  layer).  The  summer 
subsurface  mean  temperature  fields  arc  quite  different  from 
the  mean  SST  (multilayer).  Furthermore,  the  summer  mean 
SST  field  is  quite  uniform  horizontally  in  the  Yellow  Sea  shell, 
whereas  the  subsurface  temperature  fields  show  a  cold  core  in 
the  central  bottom  of  the  Yellow  Sea.  i.c..  YSCW.  In  winter  it 
merges  with  the  surface  water,  but  from  April  to  November  it 


Figure  6.  Number  iff  temperature  profiles  in  the  Yellow  Sea 
shell  during  January-June  1977.  This  is  typical  data  density  in 
any  given  (>  months  during  1903-197?. 
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Figure  7.  Surface  (it  =  0)  mean  temperature  fielii  during  (a)  winter  and  (b)  summer. 


is  a  distinct  water  mass  with  temperature  of  ->  -5'C.  During 
summer  the  overlying  waters  are  10  - 15  C  warmer. 

5.  Space/Time  Sorting 

At  present  the  spatial  and  temporal  scales  are  important 
input  data  for  the  L'.S.  Navy's  Optimum  I'licrmal  interpolation 
System  (OTIS)  model:  that  is.  the  user  should  provide  the 
values  of  these  scales  before  running  the  OTIS  model  What 
are  these  values  for  the  Yellow  Sea?  We  should  use  the  most 
extensive  MOODS  data  to  obtain  these  scales.  On  the  basis  of 
the  bowl-ivpe  bathymetry  (Figure  2).  Guati  [  1 004]  suggested 
that  the  basin  may  be  treated  as  an  isolated  system.  Therefore 
it  is  reasonable  to  assume  that  the  temperature  ACFs  inside 
the  Yellow  Sea  are  isotropic;  that  is.  ACFs  only  depend  on  the 
distance  between  two  temperature  locations.  The  major  reason 
lor  using  horizontal  ispirbpicity  here  is  the  data  sparse  ness 
(Tabic  I  ).  Il  we  do  not  use  this  assumption,  the  number  of  bins 
becomes  very  large;  for  example,  it  will  be  27.000  if  each  of  the 


(a)  Winter 


Figure  X.  Middepth  (<r  =  0.5 >  mean  temper 


temporal  and  spatial  (  v  and  v)  lags  has  30  bins.  It  is  difficult  to 
obtain  meaningful  statistics. 

For  each  observation  at  certain  a  levels  ur  =  0.  0,5.  ij.s. ). 
r,;’\  find  the  closest  and  point  climatological  value  j) " 1  and 
compute  the  anomaly  T'/,’" .  by  subtracting  P/r>  from  Ty/\ 
i.e..  7s/, M  r/r  -  7;'".  Every  individual  anomaly.  7'/”'.  is 
paired  with  the  other  data  points,  TYl 1  .  within  the  four  sea¬ 
sons.  The  temporal  and  spatial  differences  or  lags  are  calcu¬ 
lated  between  the  two  anomalies.  The  anomaly  pair  (  T\“ 

' )  is  then  placed  in  the  corresponding  temporal  (with 
increment  A/  -  I  day)  and  spatial  (with  increment  A r  =  It) 
km)  lag  bin.  If  the  lags  arc  within  Ar,,  (5  km)  and  At,,  (0.5  day), 
the  corresponding  pair  is  placed  into  bin  (0. 0).  If  the  horizon¬ 
tal  lag  is  between  mAr  -  Ar,,  and  mSr  +  Ar()  and  the 
temporal  lag  is  between  //Ar  A/ ,,  and  n  St  -  Ar0l  the  pair 
is  placed  into  the  bin  (//;.  /;).  Tile  pair-number  distributions. 
Pint.  it),  for  the  four  seasons  are  depicted  in  Figure  10.  We 
see  some  uneven  distribution  in  the  temporal  and  spatial  bins. 


(b>  Summer 


ature  field  during  la)  winter  and  (b)  summer. 
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(a)  Winter 


{ b)  Summer 


Longitude  (B)  Longitude  (E) 

Finn  re  9.  Near-bottom  (cr  =  O.S)  mean  temperature  field  during  (a)  winter  and  (b)  summer. 


Almost  everywhere.  /'( m .  n )  is  equal  to  or  more  than  51  It).  The 
maximum  /*(>».  w)  is  located  in  bins  with  near  1-2  day  tem¬ 
poral  l.'tgs  and  SO— 1 20  km  spatial  lags. 
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6.  Autocorrelation  functions  which  varies  with  the  spatial  and  temporal  lags  (m.  n)  and 

6.1.  ACF  Estimation  level  tr.  Surface  ACFs.  ij<n,(  w. //).  of  both  winter  (Figure  1  la) 

After  the  anomaly  pairs  have  been  spatially  and  temporally  and  summer  (Figure  I  lb)  are  characterized  by  decreasing  am- 

sorted.  the  ACF  value  for  each  bin  (in.  it)  is  computed  by  plilude  with  longer  temporal  and  spatial  lags. 
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Figure  10.  The  pair-number  (in  100)  distribution  in  the  (m,  n )  space  for  (a)  spring,  (b)  summer,  (c)  fall,  and 
(d)  winter. 
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6.2.  The  t  Test  for  the  AC'F  Estimation 

Equation  (3)  indicates  that  the  computed  ACE  in  the  bin 
lm.nl  is  in  fact  the  estimation  of  the  correlation  coefficient  of 
pairs  in  that  bin.  A  t  value  for  verifying  the  significance  of  a 
sample  correlation  coefficient  is  given  bv  \  Walpole  and  Sixer s. 
1989] 


which  is  a  value  of  the  statistic  T  having  a  t  distribution  with 
P  -  2  degree  of  freedom.  Using  a  as  the  level  of  significance, 
a  criterion  t)„  is  obtained 


The  criterion  p,  has  a  strong  seasonal  variation  (Figure  12) 
due  to  the  change  of  pair  numbers  P(ni.  n).  The  criterion  tjm 
is  much  smaller  for  most  bins  in  summer  than  in  winter. 

When  i)tw>(w.  ")  >  V,.-  the  estimated  ACF  is  significant 
on  the  level  of  o.  Since  both  n‘ "’(>«,  n )  and  ij„  have  seasonal 
variations,  the  significance  of  the  ACF  estimation  should  also 
change  with  seasons.  The  significant  surface  ACF  estimation 
(o  1  0.10)  is  limited  to  the  left  lower  corner  of  the  (-.  r)  plane 
with  r  <  15  day  and  r  <  201)  km  in  winter  (Figure  13a)  and 
occupies  nearly  the  whole  left  part  of  the  ( r.  r)  plane  with  r  < 
20  day  in  summer  (Figure  13b).  Such  a  seasonal  variability  in 
significant  ACF  estimations  leads  to  a  significant  difference  in 
statistical  parameters  such  as  decorrelation  scales. 


Figure  II.  Dependence  of  surface  ACF  on  spatial  and  tem¬ 
poral  lags  On.  n)  for  (a)  winter  and  (b)  summer. 


Figure  12.  Criterion  (0.10  significant  level)  for  surface  ACF 
estimation:  (a)  winter  and  (b)  summer. 


63.  Signal-to-Noise  Ratio 

The  measured  variance  s'  of  the  thermal  fields  is  separated 
into  signal  and  noise  whereby 

sz  =  if  4-  i* 

The  noise  variance  is  brought  on  from  two  sources,  geophysical 
and  instrumentation  errors.  Here  the  geophysical  error  is  un¬ 
resolved  thermal  variability  with  scales  smaller  than  the  typical 
timescales  and  space  scales  between  two  temperature  profiles. 
In  this  study  the  unresolved  scales  arc  0.5  day  and  5  km.  The 
ACF  value  at  the  first  bin  (0.  0)  does  not  represent  the  corre¬ 
lation  between  profiles  paired  by  themselves  and  therefore 
docs  nol  equal  I.  Following  SprinJaU  and  Shyers  [1991],  the 
signal-to-noise  ratio  is  computed  by 

_  .v.  Tj(0. 0) 

ASs,  “  \l  -  Tj<0.  0)  (;>) 

The  larger  the  A  is,  the  less  geophysical  error  existed.  If 
rj(0.  0)  =  1,  there  is  no  noise.  A  =  *,  and  if  p(0.  0)  =  0,  there 
is  no  signal.  If  A  >  2.  the  ratio  of  the  signal  variance.  s~,  to  the 
noise  variance,  s',  is  greater  than  4.  which  was  considered 
quite  good  conditions  by  White  el  al.  1 1982]  and  Sprinndl  and 
Meyers  [1991]. 

6.4.  Temporal  Dependence  of  ACFs 

The  temporal  dependence  of  ACF.  V  "'(rn.  n ).  can  be 
easily  discussed  by  tj-«  curves  at  several  different  spatial  lags, 
e.g..  m  -  0  (“no  lag  "),  m  =  1  (10-krn  lag),  and  m  —  15 
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Figure  13.  Significant  surface  ACF  estimation  on  level  of 
0.10  for  (a)  winter  and  (b)  summer. 


(150-km  lag).  These  curves  were  plotted  for  different  seasons 
and  three  depths  (cr  -  U.  0.5.  0.8)  in  order  to  see  the  seasonal 
and  vertical  variations.  The  three  winter  curves.  tj(o,|(0,  n). 
VT>(  1 .  ii ),  and  tj ( 15,  n),  are  plotted  for  the  surface.  <r  =  0 
(Figure  14a).  the  middepth.  <r  -  0.5  (Figure  14b),  and  the  near 
bottom,  it  =  0.8  (Figure  14c).  We  notice  from  Figure  12a  that 
tj,0,(15.  //)  is  significant  only  if  the  temporal  lag  is  smaller 
than  10  days  (i.e..  n  <  10).  The  three  summer  curves. 

n),  ft),  and  15,  //).  arc  plotted  for  the 

surface,  o-  =  0  (Figure  15a).  the  middepth,  it  =  0.5  (Figure 
15b).  and  the  near  bottom,  cr  =  0.8  (Figure  15c). 

The  ACFs  have  the  following  features: 

1.  Its  temporal  variability  weakens  as  the  spatial  lag  in¬ 
creases  and  becomes  extremely  small  (near  constantly  low 
ACF)  at  the  spatial  lag  m  =  15  (150  km)  except  for  the 
summer  surface  field  (Figure  15a).  where  ACF  decreases  qua- 
si-linearly  with  the  lime  lag  (doited  line  in  Figure  15a). 

2.  Its  vertical  variability  is  quite  small  during  winter 
(Figures  14a.  14b.  and  14c)  and  not  so  small  during  summer 
(Figures  15a.  15b.  and  15c).  This  coincides  with  the  single-layer 
structure  in  winter  and  the  multilayer  structure  in  summer 
(P.  C.  Chu  et  al..  A  feature  model  for  the  Yellow  Sea  shelf 
thermal  structure,  submitted  to  Journal  of  Geophysical  Re¬ 
search.  10%)  (hereinafter  referred  to  as  Chu  et  al..  submitted 
manuscript.  1996). 

3.  During  summer  the  ACF s  temporal  variability  weakens 
with  depth.  The  surface  ACF  shows  a  fast  reduction  with  time 
lag  n  (Figure  15a).  The  midvvater  (<r  =  0.5)  ACF  has  slower 
reduction  than  the  surface  as  n  increases  (Figure  15b).  The 


near-bottom  (cr  =  0.8)  ACF  fluctuates  around  certain  values 
(0.5  for  no  spatial  lag.  0.28  for  10-km  lag.  and  almost  0  for  150 
km)  as  n  increases  (Figure  15c). 

6.5.  Spatial  Dependence  of  ACFs 

The  spatial  dependence  of  ACF.  7 n ).  can  lie  easily 
discussed  by  r\-m  curves  at  several  different  temporal  lags,  e.g.. 
n  =  0  ("no  lag"),  n  =  1  ( 1  -day  lag),  and  n  =  15  (15-day  lag). 
These  curves  were  plotted  for  different  seasons  and  three 
depths  (it  =  0. 0.5. 0.8)  in  order  to  sec  the  seasonal  and  vertical 
variations.  The  three  winter  curves.  0).  tj <u,(m.  I ). 

and  15).  are  plotted  for  the  surface.  <r  =  0  (Figure 

16a),  the  middepth,  a  -  0.5  (Figure  16b).  and  the  near  bottom. 
it  =  0.8  (Figure  1 6c).  We  notice  from  Figure  I2a  that  tj(0)(ot,  15) 
is  significant  only  if  the  spatial  lag  is  smaller  than  50  km  (i.e.. 
m  <  5).  The  three  summer  curves,  tj<1,,(/m.  0),  ri"n(rii.  1), 
and  rj'"’(/7i,  15).  are  plotted  for  the  surface,  it  =  l)  (Figure 
1 7a),  the  middepth,  it  =  0.5  ( Figure  1 7b).  and  the  near  bottom. 
it  =  0.8  ( Figure  1 7c). 

The  ACFs  have  the  following  features: 

1.  Its  spatial  variability  weakens  as  the  temporal  lag  in¬ 
creases  and  becomes  very  small  (smaller  than  0.2)  at  the  time 
lag  n  =  15  dav  except  for  the  summer  near-bottom  field 
(Figure  17c).  where  the  ACF’s  horizontal  variability  at  n  =  15 
days  is  quite  close  to  that  at  n  ~  0  (no  time  lag)  and  n  =  1  day 
lag.  as  shown  in  Figure  17c.  This  indicates  that  during  summer 
the  tidal  effects  (on  the  timescale  less  or  equal  than  1  day)  arc- 
important  for  the  Yellow  Sea  bottom  thermal  field. 

2.  Its  vertical  variability  is  quite  small  during  winter  (Fig¬ 
ures  16a.  16b,  and  16c)  and  not  so  small  during  summer 
(Figures  I7a.  17b,  and  17c).  This  coincides  with  the  single-layer 
structure  in  winter  and  the  multilayer  structure  in  summer 
(C'hu  el  al..  submitted  manuscript.  1996). 

3.  During  summer  the  ACF s  spatial  variability  strengthens 
with  depth.  The  surface  ACF  shows  a  relatively  weak  reduction 
versus  spatial  lag  (Figure  17a).  The  midwater  (it  =  0.5)  ACF 
has  a  stronger  reduction  than  the  surface  as  the  spatial  lag 
increases  (Figure  17b).  The  near-bottom  (tr  =  0.8)  ACF  has 
the  strongest  reduction  versus  spatial  lag  (Figure  17c). 


7.  Temporal  and  Spatial  Decorrelation  Scales 

7.1.  The  Gaussian  Model 

The  U.S.  Navy's  OTIS  uses  the  Gaussian  function  to  fit  the 
ACF  [Phoebus.  1088:  Clarity.  1983] 

r)  r'(m,  n )  =  t)'” (0.  (I)  exp  [-.t;.(»iArr'  -  C'( /iAr):]  (6) 

where  p'"'(rii.  n )  denotes  the  Gaussian  tit  of  the  ACF  value 
(at  level  it)  in  the  bin  with  the  spatial  separation  mlr  and  the 
temporal  separation  u At.  A r  and  At  are  increments  for  the 
space/time  separation.  A„  1  and  C„  1  are  horizontal  and  tem¬ 
poral  decorrelation  scales  at  different  levels  it.  The  dominant 
space  scales  and  timescales  obtained  from  17  are  important  not 
only  in  determining  sampling  density  but  also  for  the  optimum 
interpolation  of  the  observed  data.  The  U.S.  Navy's  OTIS  runs 
everyday  and  combines  real-time  data,  climatology,  and  pre¬ 
dictions  from  ocean  mixed  layer  models  to  represent  an  accu¬ 
rate  picture  of  the  ocean  thermal  structure  on  global  and 
regional  scales  [Phoebus.  1988:  Clancy.  1983].  Before  running 
OTIS,  the  user  should  specify  temporal  and  spatial  scales .4., ' 
and  C., 1 .  whose  values  are  well  determined  in  deep  waters  but 
not  in  continental  shelves.  Therefore  determination  of  tempo- 
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Figure  14.  Winter  temporal  variation  of  ACF  at  different  spatial  lags:  m  =  (I  (no  lag.  dash-dotted  line),  nt  ~ 
1  (10-km  lag.  solid  line),  and  m  =  15  ( 150  km,  dotted  line)  for  three  levels:  (a)  surface  (a-  =  0).  (b)  middepth 
(rr  =  OJ),  and  (c)  near  bottom  (rr  =  0.8). 


ral  and  spatial  decorrelation  scales  is  crucial  for  the  continental 
shelf  thermal  interpolation. 

7.2.  The  F  Test  for  the  Gaussian  Model 
One  can  test  the  hypothesis  H„  that  the  Gaussian  model  (6) 
is  not  significant  by  merely  forming  the  ratio  [ Walpole  and 
Myers,  1989] 

SSRfk 

f~SSE/{l  -  k  -  1)  (7) 

where  k  -2.1  is  the  total  number  of  bins  in  either  spatial  or 
temporal  lags,  and 

SSR  =  £  X  [In  fT'fm.  n)  -  In  V"  ]\ 

m  n 

(«) 

SSE  =  X  [1°  n)  -  In  tj(/w.  «)]: 

«  n 

denoting  a  regression  sum  of  squares  and  a  residual  mean 
square,  respectively.  When 


/>/.(*,/“*  -  1> 

we  reject  //„.  Here /„(/:,  /  -  k  -  1)  satisfies  the  F  distribu¬ 
tion  with  v,  =  k,  —  l  —  k  —  1,  and  a  confidence  level  of 
a.  The  f  values  for  different  levels  and  seasons  are  listed  in 
Table  2.  All  the  values  of  /  exceed  the  critical  value  (5.45)  of 
the  F  distribution  for  2  and  28  degrees  of  freedom  at  a  -  0.01. 
Therefore  we  can  confirm  that  the  Gaussian  model  is  reason¬ 
able  for  the  Yellow  Sea  thermal  ACF. 

7.3.  Seasonal  Variabilities  of  the  Decorrelation  Seales 

The  computed  ACFs  for  different  seasons  and  levels  (it  =  U. 
0.5,  0.8)  were  then  fitted  to  a  Gaussian  function  of  the  form  of 
(6)  by  the  regression  method,  which  leads  to  the  spatial  and 
temporal  decorrelation  scales.  .4.) 1  and  C„  \  respectively.  Ta¬ 
ble  3  presents  these  decorrelation  scales  and  sigr»al-to-noise 
ratios  for  winter  and  summer. 

The  signal-to-noise  ratios  for  both  winter  and  summer  are 
higher  than  that  for  the  deep  waters  of  the  eastern  equatorial 
Pacific,  which  is  around  1.0  as  reported  by  Sprintall  and  Meyers 
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Figure  15.  Summer  temporal  variation  of  ACF  at  different  spatial  lags:  m  —  0  (no  lag.  dash-dotted  line). 
m  =  1  (10-km  lag,  solid  line),  and  m  =  15  (150  km.  dotted  line)  for  three  levels:  (a)  surface  (r r  -  0),  (b) 
middepth  (<r  =  0.5).  and  (e)  near  bottom  (a-  =  0.8). 


(1991],  This  infers  that  the  temperature  signal  in  the  Yellow 
Sea  shelf  is  stronger  than  in  the  open  waters,  in  both  winter 
and  summer  seasons  the  largest  noise  occurs  at  the  midlevel 
(«•  =  0.5).  This  is  expected  since  this  is  the  transition  layer 
where  both  the  meteorological  and  topographic  effects  occur. 
The  signal-to-noisc  ratio  is  greater  in  summer  than  in  winter. 

The  vertical  variation  in  temporal  and  spatial  decorrelation 
scales  is  smaller  in  winter  than  in  summer.  In  winter  the  tem¬ 
poral  scale  varies  only  a  half  a  day  and  the  horizontal  scale 
changes  only  14  km  among  three  different  levels  (<r  =  0.  0.5. 
0.8).  This  vertically  quasi-uniformity  in  decorrelation  scales 
also  represents  the  w  inter  single-layer  structure  caused  by  the 
strong  surface  forcing.  In  winter  the  Yellow  Sea  shelf  experi¬ 
ences  strong  mixing  due  to  both  strong  winds  and  the  upward 
buoyancy  flux.  With  a  large  part  of  the  Yellow  Sea  having 
depths  less  than  5(1  m.  vertical  mixing  reaches  the  bottom  and 
creates  isothermal  profiles.  Thus  decorrelation  scales  will  be 
similar  front  the  surface  to  the  near-bottom  (<r  =  0.8)  water 
column.  In  summer  the  temporal  scale  increases  5  days  and  the 


horizontal  scale  decreases  04  km  from  the  surface  (rr  =  0)  to 
the  near-bottom  (<r  -  0.8)  waters.  This  vertically  varying  de- 
correlation  scale  also  implies  the  summer  multilayer  structure. 

Surface  horizontal  decorrelation  settles  arc  almost  100  km 
longer  in  summer  than  in  winter.  This  is  due  to  the  strong  solar 
heating  in  summer,  causing  a  relatively  uniform  SST  field. 
Surface  temporal  decorrelation  scales  are  2.4  days  shorter  in 
summer  than  in  winter.  This  might  be  caused  by  the  shallower 
surface  mixed  layer  in  summer  (less  thermal  inertia).  Only  the 
upper  layer  water  is  affected  in  summer  by  the  atmospheric 
forcing  rather  than  the  entire  water  column  as  in  winter.  The 
surface  temporal  decorrelation  scale  should  be  shorter  in  sum¬ 
mer  than  in  winter. 

An  interesting  feature  shown  in  Table  3  is  the  increase  of 
temporal  decorrelation  scale  with  depth  in  both  summer  (ev¬ 
ident)  and  winter  (slight).  The  near-bottom  water  (<t  -  0.S)  has 
the  longest  temporal  scale  in  summer,  which  could  be  directly 
related  to  the  existence  of  the  Yellow  Sea  Cold  Water  through¬ 
out  the  summer  in  the  middle  of  the  Yellow  Sea. 
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Figure  16.  Winter  spatial  variation  of  ACT  at  different  temporal  lags:/i  =  0  (no  lag.  dash-dotted  line),  n  - 
I  (1-day  lag.  solid  line),  and  ;;  =  15  (15-day  lag.  dotted  line)  for  three  levels:  (a)  surface  (;r  =  0).  (b)  middepth 
(<r  =  0.5).  and  (e)  near  bottom  (a  =  0.8). 


8.  Applications 

8.1.  Optimum  Interpolation 

Optimum  interpolation  assigns  a  weight  to  each  observation 
that  accounts  for  variation  in  spatial  and  temporal  sampling. 
The  interpolated  temperature  anomaly  at  the  grid  point.  T',-.  is 
a  linear  combination  of  the  observed  anomalies.  T] .  T\.  .  T\ 
with  weights  a,,  cr-.  ••• .  «* 

S 

f;,=  S«  J’  (y) 

t-i 

The  values  of  the  weights.  u,(;  =  1.2.***.  :V).  are  found  by 
minimizing,  in  a  least  squares  sense,  the  difference  between  the 
interpolated  value  at  the  grid  node  and  the  tme  value  there  and 
obtained  from  solving  the  algebraic  equations  \Gamlin.  into]: 

•v 

^  «,p,„  +  A  or,  =  n,7l  (10) 

n-l 


where  ;  =  l,  2,  -  •  • ,  ,V.  Here  A  is  the  signal-to-noise  ratio. 
is  the  autocorrelation  between  locations;  and  j.  and  's  ■he 
autocorrelation  between  the  grid  node  and  location  /. 

=  »] l m„,  «„).  Mn,  -  n,„)  (11) 

where  (»>„.  n tj )  and  »a,i  represent  spatial  and  tem¬ 

poral  separations  between  locations;  and  j  and  between  k  na¬ 
tion  i  and  the  grid  node,  respectively.  For  each  grid  node 
location.  ( 10)  results  in  a  set  of  .V  linear  algebraic  equations  to 
be  solved  for  V  unknowns.  by  matrix  inversion.  In  ( 10)  the 
autocorrelations  p,(J  and  n(il  and  the  signal-to-noise  ratio  A 
have  been  computed  at  three  levels  (<r  -  0.  0.5,  0.8)  from  the 
MOODS  data  set.  Therefore  for  these  levels  the  V  weights  a, 
can  be  determined  by  solving  (10),  and  therefore  any  new 
observations  T\.  T'..  ••*.  T\  can  lx*  interpolated  into  any 
grid  node  G. 

To  run  the  U.S.  Navy's  OTIS,  the  temporal  and  spatial 
decorrelation  scales  .-I  „  1  and  C~  1  should  be  prescribed  [Phoe- 
hu\.  1988:  Clancy .  1983].  Table  3  provides  useful  information 
for  assigning  these  values  for  the  Yellow*  Sea  OTIS  model. 
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Disurnc*  {km) 


Figure  17.  Summer  spatial  variation  of  ACT  at  different  temporal  lags:  n  —  0  (no  lag.  dash-dotted  line). 
n  =  I  { 1  -day  lag.  solid  line),  arid  «  =  15  ( 15-day  lag.  dotted  line)  for  three  levels:  (a)  surface  (w  =  0),  (b) 
middepth  (it  -  (1.5),  and  (c)  near  bottom  (<r  =  0.8). 


8.2.  Observational  Network  Design 

The  temporal  and  horizontal  decorrelation  scales  are  useful 
for  an  optimal  observational  network  design.  As  we  mentioned 
in  section  6.3.  the  noise  comes  from  instrumental  and  geo 
physical  errors.  Since  the  instrumental  error  in  expendable 
bathythermograph  (XBT)  measurements  is  usually  about  0.1 'C' 
[Harnett  and  Bernstein,  1980)  and  Is  even  smaller  in  conductiv- 
ity-ternperature-depth  (CTD)  measurements,  the  instrumental 
error  is  generally  neglected  against  the  geophysical  error.  This 
implies  that  the  curtailment  of  noise  must  be  accomplished  by 
the  reduction  of  geophysical  error.  This  is  usually  done  by 
increasing  the  sample  density.  Having  determined  the  statisti¬ 
cal  structure  of  thermal  variability  in  the  Yellow  Sea  shell,  the 


Table  2.  The/  Values  of  the  Yellow  Sea  ACF  Intimations 

Winter  Summer 

Surface  28.64  12.26 

0.5  26.43  11.71 

0.8  15.01  11.34 


minimum  sampling  density  required  to  detect  thermal  variabil¬ 
ity  can  now  be  arbitrated  as  two  or  three  samples  per  decor¬ 
relation  scale  \Spnntall  and  Meyers.  1991  j.  This  would  mean 
that  spatially,  any  temperature  measurement  in  both  summer 
and  winter  may  he  conducted  at  50-S0  km  and  4-6  day  inter¬ 
vals  with  the  knowledge  that  the  subsurface  features  will  also 
be  adequately  sampled. 


Table  3.  Seasonal  and  Vertical  Variations  of  Decorrelation 
Scales 


Season 

Depth 

Temporal 

Decorrelation 

Scale. 

days 

Spatial 

Decorrelation 

Scale. 

km 

Signal -fo- 
N'oise 
Ratio 
(A) 

Winter 

surface 

14.7 

158 

2.06 

50% 

14.7 

167 

1.88 

Stffi 

15.2 

172 

121 

Summer 

surface 

12.3 

15] 

3.(1 

505? 

15. K 

169 

121 

sir,. 

17.2 

157 

5.0 

(Ill  HI  A I ..  SCALES  OF  THE  'i  ELLOW  SFA  THERMAL  VARIABILITY- 
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9.  Conclusions 

1.  The  computation  of  ACT  for  <r  levels  proposed  in  this 
paper  demonstrates  a  good  capability  to  obtain  the  statistical 
structure  of  the  continental  shelf  thermal  field.  After  the  sta¬ 
tistical  structures  for  several  a  levels  were  constructed  front 
35.658  profiles  ( 1929- 1991)  on  the  Yellow  Sea  shelf  from  the 
MOODS  database,  we  obtain  the  temporal  and  spatial  decor¬ 
relation  scales  and  their  seasonal  and  vertical  variations  for  the 
Yellow  Sea  shelf  thermal  fields.  These  are  important  for  opti¬ 
mum  interpolation;  for  example,  the  decorrelation  scales  arc 
crucial  for  running  the  U.S.  Navy's  OTIS  model. 

2.  The  signal-to-noisc  ratios  for  both  winter  and  summer 
are  quite  high.  This  infers  that  the  temperature  signal  in  the 
Yellow  Sea  shelf  is  strong.  In  both  the  winter  and  summer 
seasons  the  largest  noise  occurs  at  the  midlevel  (<r  =  (l_5).  The 
signal-to-noisc  ratio  is  greater  in  summer  than  in  winter. 

3.  The  vertical  variation  in  temporal  and  spatial  decorre¬ 
lation  scales  ts  smaller  in  winter  than  in  summer.  In  winter  the 
temporal  scale  varies  only  a  half  a  day  and  the  horizontal  scale 
changes  only  14  km  among  three  different  levels  (tr  =  0.  0.5, 
0.8).  This  vertical  quasi-uniformity  in  decorrelation  scales  also 
represents  the  winter  single-layer  structure  caused  by  the 
strong  surface  cooling  and  stirring.  In  summer  the  temporal 
scale  increases  5  days  and  the  horizontal  scale  decreases 94  km 
from  the  surface  (tr  =  0)  to  the  near-bottom  (<r  =  0.8)  waters. 
This  vertically  varying  decorrelation  scale  also  implies  the  sum¬ 
mer  multilayer  structure. 

4.  Tlte  seasonal  variation  in  temporal  and  spatial  decorre¬ 
lation  scales  has  its  maximum  at  the  surface,  reduces  with 
depth  to  a  minimum  at  midlevel  (<r  =  0.5).  and  then  enhances 
with  depth  to  the  near-bottom  level.  The  surface  horizontal 
decorrelation  scales  are  almost  100  km  longer  in  summer  than 
in  winter.  This  is  due  to  the  strong  solar  heating  in  summer, 
causing  a  relatively  uniform  SST  lield.  Surface  temporal  deco¬ 
rrelation  scales  are  2.4  days  shorter  in  summer  than  in  winter. 
This  might  be  caused  by  the  shallower  surface  mixed  layer  in 
summer  (less  thermal  inertia). 

5.  The  temporal  decorrelation  scale  increases  with  depth 

in  both  summer  (evident)  and  winter  (slight).  The  near-bottom 
water  (<r  0.8)  has  the  longest  temporal  scale  in  summer, 

which  could  he  directly  related  to  the  existence  of  the  Yellow 
Sea  Cold  Water  (YSCVV)  throughout  the  summer  in  the  mid¬ 
dle  of  the  Yellow  Sea. 

6.  Tlte  temporal  and  horizontal  decorrelation  scales  are 
useful  for  designing  an  optimum  observational  network.  Tlte 
minimum  sampling  density  required  to  detect  thermal  variabil¬ 
ity  in  the  Yellow  Sea  shell  would  he  50-  80  km  and  4-6  day 
intervals  per  temperature  measurement  with  the  knowledge 
that  the  subsurface  features  w  ill  also  be  adequately  sampled. 
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